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Abstract. We present observations of 11 quasars, selected in the range z m 2.2-4.1, obtained with ESA's 
Superconducting Tunnel Junction (STJ) camera on the WHT. Using a single template QSO spectrum, we show 
that we can determine the redshifts of these objects to about 1%. A follow-up spectroscopic observation of one 
QSO for which our best-fit redshift (z = 2.976) differs significantly from the tentative literature value [z ~ 2.30) 
confirms that the latter was incorrect. 
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1. Introduction 

Large ground and space telescopes combined with solid 
state detectors have revolutionized optical astronomy over 
the past two decades, yet deriving physical diagnostics 
of stars and galaxies still requires the somewhat indirect 
methods of filter photometry or dispersive spectroscopy to 
measure spectral features, energy distributions, and red- 
shifts. The recent development of high-efficiency supercon- 
ducting detectors ( Perryman et al., 1993 ; Peacock et al. 
|1996| ) has introduced the possibility of measuring indi- 
vidual optical photon energies directly, and the first high 
time-resolution spectrally-resolved observations of rapidly 
variable sources such as cataclysmic variables and op- 
tical pulsars using these techniques have been reported 

Romani et al 



1999 



1999; Perryman 



(Perryman et al. 
[et al., 2001 ; Bridge et al., 2001). Many extensive obser- 



vational programmes which aim at determining the large- 
scale structure of the Universe, and galaxy formation and 
evolution (e.g., the Sloan Digital Sky Survey, Fan et al 



1999; the Anglo- Australian Telescope 2dF survey, Croom 



et al. 2001), demand high-efficiency extragalactic spec- 



troscopy. Here we report the first optical measurements of 
spectral energy distributions of quasars using an imaging 
detector with intrinsic energy resolution, and show that 
we can determine their redshifts directly with excellent 
precision. 
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2. Observations 

We observed 11 quasars in the redshift range z = 
2.2-4.1, the sample comprising relatively bright high- 
redshift Lyman-limit quasars from the published litera- 
ture ( Sargent et al., 1989| ), supplemented by three lower 
redshift objects, two of which were discovered in objec- 
tive prism- type surveys (Table [l]). Observations used the 
ESA superconducting tunnel junction (STJ) camera, S- 
Cam2 ( |Rando et al, 2000| ), on the 4.2-m William Herschel 
Telescope, La Palma, between 2000 October 1-4. The 
camera is a 6 x 6 array of 25 x 25 /im 2 (0.6 x 0.6 arcsec 2 ) 
tantalum junctions, providing individual photon arrival 
time accuracies to about 5 /is, a resolving power of 1Z ~ 8 
at A = 500 nm, and high sensitivity from 310 nm (the 
atmospheric cutoff) to about 720 nm (currently set by 
long-wavelength filters to reduce the thermal noise pho- 
tons). All objects show strong Ly-a and CIV emission lines 
which, at these redshifts, will be present within our wave- 
length response range. Observations were made in modest 
seeing (1-1.5 arcsec at airmass X = 1), and at air-masses 
between X = 1.07-1.82. 

3. Data reduction 

Information on each detected photon consists of arrival 
time, x, y coordinates of the junction, and an energy chan- 
nel in the range 0-255. A photon of energy E p (in eV) 
incident on a particular junction is assigned to an energy 
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Table 1. The 11 quasars observed. V gives the catalogue 
magnitude (probably questionable for 2143—158). T gives 
the exposure time in seconds, and z b s our estimated red- 
shift. The final two columns give the literature redshift 
and its s ource: S89 = ^argent et al. (198S| ) [spectroscopy] ; 
M77 = [MacAlpine et al (1977|) [objective prism]; C91 



= Chaffee et al. (1991) [spectroscopy]; C85 
et al. (1985) [grens plate] 
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channel Ei = G ■ E p + C, where each pixel is characterised 
by its own gain G (in channels per eV) and offset C (in 
channels). Laboratory measurements have confirmed that 
all 36 junctions have a highly linear, albeit slightly pixel- 
dependent, energy response. Calibration consists of first 
bringir 
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ale, corresponding to an arbitrary reference pixel, 
fixed gain map based on laboratory measurements. 

of the reference pixel is constant (C = —2.0), 
gain is then the only free parameter in the ab- 
5nergy calibration. Small temporal gain variations 
g from bias voltage drifts and small detector tem- 
'e variations (±0.01 K on the nominal operating 
ature of «0.32 K) are monitored and calibrated, 
traction of the appropriate sky contribution for 
q basar spectrum can in principle be based on the 



background signal in the outer array junctions, but given 
the small array size and seeing and refraction effects, we 
generally also took a nearby sky frame immediately follow- 
ing each quasar observation. Most observations were taken 
in astronomically dark time; QSO 2233+136, 2143-158, 
and 0148—097 were observed with the Moon setting, with 
background subtraction slightly less accurate. 

4. Results 

We have determined each quasar redshift by comparing 
the calibrated energy distributions, f bs(Ei), with a sin- 



gle rest- frame composite quasar spectrum (Zheng et al 



1997) based on 284 Hubble Space Telescope Faint Object 
Spectrograph spectra of 101 quasars with z > 0.33. For 
a given gain G and redshift z, we construct the model 
energy-channel distribution fmod{Ei), as follows. The tem- 
plate spectrum is shifted from the rest frame to redshift 
z, and a mean accumulated absorption of the Lyman for- 
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Hesults for QSO 0127+059, 0148-097, and 
Left: the observed (black curves) and mod- 
;y curves) energy channel distributions (arbi- 
ts). Insets indicate the estimated Poisson noise, 
above the top left panel show the mapping be- 
y channel and wavelength. Right: the cor- 
dependence of % 2 on z. Vertical dashed 
the literature redshifts; the dotted line for 
!7+059 indicates the spectroscopic redshift re- 
this letter [z = 3.04). 
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1990) (all our objects are at high Galactic latitude, and 
we neglect Galactic reddening) . The resulting spectrum is 
corrected for the mean atmospheric transmission at the 
relevant airmass, adjusted for the instrument and tele- 
scope efficiency curves and exposure time, transformed 
from wavelength spectra to energy-channel spectra, and 
finally convolved with a suitable Gaussian in order to ac- 
count for the finite energy resolution of the detector. We 
then derive redshift z and gain G (and a normalization 
constant N), by minimizing the classical \ 2 function: 



170 r 



i=50 



fobs(Ei) — N ■ f mo d(Ei) 



n 2 



a f<,UEi) 



(1) 



est for this redshift is introduced (M0ller & Jakobsen, 



using a downhill simplex routine (Press et al., 1995). 
Summation extends over relevant energy channels, and 
118 is the number of degrees of freedom. Resulting gains 
are in the range 42-45 channels eV -1 , consistent with lab- 
oratory calibration. The resulting redshifts are listed in 
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Fig. 2. Observed versus literature redshifts. Numbers re- 
fer to the objects listed in Table |], and symbol sizes cor- 
respond to x 2 (smaller symbols indicating a poorer fit). 
QSO 0127+059 has an incorrect literature redshift of 2.30; 
our spectroscopic follow-up observation yields z = 3.04, 
moving the point to the position shown in grey. The 
dashed line shows the nominal 1:1 correlation. 

Table [l], along with the literature values. Examples of the 
observed and modeled spectra are shown in Figure 1. The 
overall shape of these spectra, and in particular the sharp 
falloff at low energy channels, is due to the response of 
the instrument and telescope. In practice, the Ly-a emis- 
sion line and the associated break at shorter wavelengths 
contribute most to the redshift determination. Figure ^| 
compares the best-fit redshifts with the literature values. 

QSO 0127+059 is our single prominent outlier. It was 
discovered in a thin prism survey (MacAlpine et al., 1977), 
classified as a possible quasar, and tentatively assigned 
a redshift of z ss 2.30, but with an uncertain line iden- 
tification. Although the quality of our fit is acceptable 
(Figure 1), our derived redshift, z = 2.976, differs signifi- 
cantly from the literature value. We subsequently obtained 
a 1200 s spectrum of QSO 0127+059 (Figure §) with the 
Siding Spring Observatory 2.3-m telescope. The wave- 
length coverage (not optimised for quasar spectroscopy) 
was 345-537 and 560-753 nm, using the Double Beam 
Spectrograph with dichroic 1 and gratings 600B and 600R. 
We determine a spectroscopic redshift z = 3.04, which 
agrees with our estimate to about 2% (Figure ||). 

A small systematic offset in the overall correlation, of 
«0.03 in z, can be attributed to a small mismatch in the 
shape of the energy broadening function or the overall 
throughput used in the modeling. The mean scatter for 
all observations is a z — 0.03; removing the systematic 
offset, 8 of the 11 objects agree to within 1%. Several fac- 
tors, such as gain variations, erroneous sky subtraction 
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Fig. 3. The spectrum of QSO 0127+059 obtained with the 
Siding Spring Observatory 2.3-m telescope, and smoothed 
with a 15 A FWHM Gaussian. We determine z = 3.04; 
the resulting redshifted locations of Ly-a (121.6 nm), NV 
(124.0 nm), and CIV (154.9 nm) are indicated. 

or extinction correction (e.g., due to unmodeled seasonal 
Saharan dust in the atmosphere), or template mismatch 
at the object level (related to continuum slope, line ratios, 
etc.), may explain the observed spread. Formally, none of 
the fits is particularly good, in the sense that none of them 
has reduced % 2 « 1. A key factor in x 2 -statistics, however, 
is the absence of systematic errors, which will exist here 
in part due to template mismatch, although details are 
largely hidden as a result of the limited detector resolu- 
tion. The general consistency between the models and the 
observations, and the pronounced, deep and narrow, min- 
ima in all x 2 versus z plots, nonetheless indicate that our 
fits, as a set, are acceptable. 

The pronounced minima are apparent in our data sets 
truncated a posteriori to observation times as small as, 
e.g. 10-20 s for the z = 4.1 object QSO 0000-263, where 
~350 source photons s _1 were recorded. 

5. Discussion 

Although extraction of detailed physical information from 
the spectra is limited by the modest resolving power 
(1Z ~ 8) of the current device, a significant improvement in 
energy resolution can be expected in the future ( Perryman 
|et al., 1993 : Peacock et al., 1997), and additional template 



spectra could then be used for model fitting. Our results 
show that low-resolution spectroscopy of faint extragalac- 
tic sources is possible with these devices, enabling the de- 
termination of redshift, and perhaps morphological type 



and emission and absorption line ratios (lakobsen, 1999 



Mazin & Brunncr, 2000). Our instrument development is 
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aimed at larger format arrays to facilitate sky subtrac- Zheng, W., Kriss, G. A., Telfer, R. C, Grimes, J. P., & 
tion and possibly for multi-object spectroscopy, and an Davidsen, A. F. 1997, ApJ, 475, 469 
increased energy resolution to improve physical diagnos- 
tic capability. An overall wavelength response extending 
further to the red, consistent with the fundamental de- 
vice response characteristics, would also open up a larger 
accessible redshift range. 
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